Laser probing of atomic Cooper pairs 
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We consider a gas of attractively interacting cold Fermionic atoms wfiidi are manipulated by 
laser ligfit. The laser induces a transition from an internal state with large negative scattering 
length to one with almost no interactions. The process can be viewed as a tunneling of atomic 
population between the superconducting and the normal states of the gas. It can be used to detect 
the BCS-ground state and to measure the superconducting order parameter. 
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Successful cooling of trapped Fermionic ''°K atoms 
down to temperatures where the Fermi degeneracy sets 
in was reported recently j^,^. This breakthrough opens 
up new opportunities for studying fundamental quantum 
statistical and many-body physics. A major advantage 
of Fermionic atoms compared to electrons in condensed 
matter is the richness of their internal energy structure 
and the possibility to accurately and efficiently manip- 
ulate these energy states by laser light and magnetic 
fields. Furthermore, atomic gases are dilute and weakly 
interacting thus offering the ideal tool for developing and 
experimentally testing theories of many-body quantum 
physics. A major goal is to observe the predicted 
BSC-transition for Fermionic atoms - this would compare 
to the experimental realization of atomic Bose-Einstein 
condensates |^. It is still, however, an open question 
how to observe the BCS-transition, because the value of 
the superconducting order parameter (gap) is expected 
to be small and electro-magnetic phenomena such as the 
Meissner effect do not take place. 

We propose a method to detect the existence of the 
BCS-ground state and to measure the gap using laser 
light. The main difference in our method compared to 
the proposals of using off-resonant light scattering [^|J^] 
is that the light is resonant, that is, population is trans- 
ferred from one internal state to another. Furthermore, 
the final state of this process is chosen to have negligible 
scattering length compared to the initial one: this effec- 
tively creates a superconducting - normal state interface 
across which the atomic population can move. There is 
a conceptual analogy to electron tunneling from a super- 
conducting metal to a normal one, although the physical 
systems and their describtions differ in an essential way. 
Also non-optical phenomena, such as collective and sin- 
gle particle excitations, have been proposed to be used 
for observing the BCS-transition The specific 

feature of our method with respect to other suggested 
probes, both optical and mechanical (modulating the 
trap frequency), is the utilization of the superconduct- 
ing - normal state interface. The advantage is that the 
normal state population can be initally zero, thus even a 



small change in it is a significant relative effect. Further- 
more, our work offers a basic starting point for describ- 
ing any phenomenon related to the superconducting- 
normal interface. For instance, our method can be un- 
derstood as an outcoupler. With small modifications, 
one could also describe outcoupling of pairs instead of 
single atoms which would create an atomic beam with 
BCS-type statistics. 




FIG. 1. Probing of the gap in a gas of attractively interact- 
ing cold Fermionic atoms, a) Laser excitation with the cou- 
pling f2 and the detuning A transfers a Cooper paired atom 
from the internal state \g) to the state |e). b) The other atom 
in the initial Cooper pair becomes an excitation in the BCS 
state, therefore the laser has to provide also the additional 
gap energy Ag. In this picture the Fermi levels /Xg and /ie for 
the internal states have been chosen to be different from each 
other but they could also equal. 

We consider atoms with three internal states available, 
say |e), \g), and \g'). They are chosen so that the inter- 
action between atoms in states \g) and \g') is relatively 
strong and attractive, and all other interactions are neg- 
ligible. The laser frequency is chosen to transfer popu- 
lation between |e) and \g), but is not in resonance with 
any transition which could move population away from 
the state \g'). If there is more or less equal amount of 
atoms in states \g) and \g') they can form the supercon- 
ducting BCS ground state; atoms in |e) are in the normal 
state. For small intensities, the laser-interaction can be 
treated as a perturbation, the unperturbed states being 
the normal and the superconducting state. The transfer 
of atoms from \g) to |e) is then analogous to tunneling 
of electrons from a normal metal to a superconductor in- 
duced by an external voltage, which can be used as a 
method to probe the gap and the density of states of 
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the superconductor [|ri| . In our case the tunnehng is be- 
tween two internal states rather than two spatial regions, 
this resembles the idea of internal Josephson-oscillations 
in two-component Bose-Einstein condensates Fig- 
ure 1 illustrates the basic idea. The observable carrying 
essential information about the superconducting state is 
the change in the population of the state |e), we call this 
the current /. Below we calculate the current first in the 
homogeneous case and then assuming that the atoms are 
confined in a harmonic trap. Finally will discuss possibil- 
ities for experimental realization of the idea, and choosing 
the states |e), \g), and \g') for ^Li. 

We define a two-component Fermion field V'(x) = 
(■06 (x) ■i/'g(x))^, where tpe and tpg fulfil standard 
Fermionic commutation relations. The fields ipe/g can be 
expanded using some basis functions (e.g. plane waves 
or trap wave functions) and corresponding creation and 

annihilation operators: Ve/g(x) = J2j '^'j^^ '^''j^^ i^)- The 
annihilation and creation operators fulfil {cf\c^} = 

and {c^^^^, c^^^} = Sij. The two components of the field, 
corresponding to the internal states |e) and \g), are cou- 
pled by a laser. This can be either a direct excitation 
or a Raman process; we denote the atomic energy level 
difference by tUa {h = 1), the laser frequency ujl and the 
wave vector fc^ - in the case of a Raman process these are 
effective quantities. In the rotating wave approximation 
the Hamiltonian reads 

H = H, + Hgg, + ^ J d3a;V't(x)a,0(x) 

Here A = uja — ujl is the (effective) detuning, and f7(x) 
contains the spatial dependence of the laser field mul- 
tiplied with the (effective) Rabi frequency. The parts 
He and Hggi contain terms which depend only on ipe or 
tpg, tjjgi , respectively. Possible spatial inhomogeneity, e.g. 
from the trap potential, is also included in and Hggi . 

The observable carrying essential information about 
the superconducting state is the rate of change in pop- 
ulation of the |e) state. We may call it, after the 
electron tunneling analogy, the current I{t) — —{Ne), 
where = J (Pxtpl{'x.)ipei^)- The current I{t) is calcu- 
lated considering the tunneling part of the Hamiltonian, 
Ht = H - (He + Hgg, + A/2{Ne - Ng)) as a perturba- 
tion; the current / becomes the first order response to 
the external perturbation caused by the laser. We cal- 
culate it both in the homogeneous case and in the case 
of harmonic confinement. The calculations are done in 
the grand canonical ensemble, therefore the chemical po- 
tentials fig and fie are introduced. Also the detuning A 
acts like a difference in chemical potentials, thus it be- 
comes useful to define an effective quantity of the form 
A = fie — fig + A = Afi -f A. In the derivation we assume 
finite temperature, but here we only quote the results for 



T = 0. 

Homogeneous case: The assumption of spatial homo- 
geneity is appropriate when the atoms are confined in 
a trap potential which changes very little compared to 
characteristic quantities of the system, such as the co- 
herence length and the size of the Cooper pairs. In the 
homogeneous case we expand the Fermion fields V'e/g i^ito 
plane waves. The Hamiltonian becomes 



H = He 
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We calculate the 



'-kl ^ yj d3xr!(x)e*-''e-*- 
current / = —{Ne) = ^i{[H, Ne]) treating Ht as a per- 
turbation: terms of higher order than H^ are neglected. 
Because we are interested in the current between the su- 
perconducting and normal states, correlations of the form 



(and h.c.) are omitted since they correspond 
to tunneling of pairs (Josephson current). The current 
can be written 



dte{t){e-^^'{[A\0),A{t)]) - e'^*([^(0), At(i)])) 



= Efc-i^fcicf (i)cf(i), and 



e""^c,'"e where K ^ H - fieNe - fJ-gNg. The two 
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terms in the above equation have the form of retarded 
and advanced Green's functions. These are evaluated us- 
ing Matsubara Green's functions techniques, which leads 
to / = Efc; \Th?Coo tMe) - np{e + A)]A,(fc,e + 
A)Ae{l,e), where np are the Fermi distribution func- 
tions and Agje are the spectral bmctions for the su- 
perconducting and normal states. We use the standard 
expressions Ae[l,e) = 27rJ(e — ^i) and Ag{k,e -I- A) = 
2Tr[uld{e+A-ujk)+vlSie+A+LJk)] Here 6 ^ Ei-fie 
and Uk , Wfc and ujk are given by the Bogoliubov transfor- 
mation. Because we are interested in the tunneling out of 
the superconductor, only the term proportional to vf. is 
considered. The laser field is chosen to be a running wave, 
that is r2(x) = Q.e'^^^'^. The term \Tki\'^ now produces a 
delta-function enforcing momentum conservation. Note 
that this is very different from the assumption of a con- 
stant transfer matrix (X]fc; l^fciP — ^ I^PX^fcz) made in 
the standard calculation for tunneling of electrons over 
a superconductor ~ normal metal surface The final 
result becomes (assuming T = 0) 
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where k is given by the following energy conservation 
condition: 



_ + = 0, Wfc - V '^fc 

and p is the density of states which appears when the 
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summation over momenta is changed into an integration 
over energies. 

The laser momentum can be very smah compared 
to the momentum of the atoms, especially in the case of 
a Raman process. By setting kL = the result becomes 
(choosing A < i.e. current from \g) to |e)) 



2A^A)^. 
^ ^A2 
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To understand the results in terms of physics, let 
us first consider the case of equal chemical potentials 
A/i = 0: 



/ = -27rl7V6'(-A - Ac 
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In order to transfer one atom from the state \g) to |e) 
the laser has to break a Cooper pair. The minimum en- 
ergy required for this is the gap energy Ag , therefore the 
current does not flow before the laser detuning provides 
this energy - this is expressed by the step function in 
(^. As |A| increases further, the current will decrease 
quadratically. This is because the case |A| = Ag corre- 
sponds to the transfer of particles with p = pp, whereas 
larger |A| means larger momenta, and there are simply 
fewer Cooper-pairs away from the Fermi surface. This 
behaviour is very different from the electron tunneling 
where the current grows as (eV)'^ — A^ (the volt- 
age eV corresponds to the detuning A in our case) be- 
cause all momentum states are coupled to each other. 

When the chemical potentials are not equal the situa- 
tion is more complicated, but the basic features are the 
same: i) treshold for the onset of the current given ba- 
sically by the gap energy, and ii) further decay of the 
current because the density of the states that can fulfil 
energy and momentum conservation decreases. 

Harmonic confinement: The trap wave functions pro- 
vide the natural basis of expansion when the atoms are 
confined in a harmonic trap. We will not, however, do 
this expansion from the beginning of the calculations but 
rather derive the current and the corresponding (position 
dependent) Green's functions directly for the Fermion 
fields tpg/g. The current becomes / = 2/m[Xret(— A)], 

t(-A) = r ^ [ d^x [ d^x'n*{^)n{^') 
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[ie(x,x',e)GL„(x',x,e + A) 
+G^,t(x,x',e- A)ig(x',x, e)] 
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and A^/g are defined Jie/g(x, x', e) — ^(G'^gf (x, x', e) — 



G^^^(x, x', e)). In order to take the imaginary part of 
the expression ^ we use the following Green's functions 



GL.(x',x,e) = 



u„(x'X(x) ^ <(xOt>„(x) 

e + Un — iS 
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Here u„(x), w„(x) and a;„ are given by the Bogohubov- 
DeGennes equations We also use the fact that the 
trap wave functions 0m(x) are real. The derivation gives 



/ = -27r^ y d^xr2(x)u„(x)(/i,„(x) 
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rf3a;rj(xX(x)(/),„(x) 



[npi-uJn) - f^F(^m)]<5(6n + A + a;„). 
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Again, we consider only the term proportional to w„, as- 
sume zero temperature and use the known form for the 
trap energy = mf2t — /ie to obtain the final result 



(i3a;rj(xX(x)0M(x) 



where the quantum number M = — A + LOn)/^t- 
Although the result is not as transparent as in the ho- 
mogeneous case, the characteristic treshold behaviour is 
expressed by the step function. Since the functions in- 
side the absolute square are typically oscillating ones, 
this term would restrict the values of n in the summa- 
tion. We have checked that with strong approximations 
one obtains again the homogeneous case result. 

Experimental realization: Typical experimental pa- 
rameters in the case of ®Li are Q.t/2n ^ IhQHz for 
the trap frequency and /ig ~ lOOkHz for the chemical 
potential [Q. The order of magnitude of the gap is 
Ac/er ~ 0.1 (see e.g. [^). In Figure 2 we show the 
dependence of the current on the detuning in the homo- 
geneous case, for Aq = 0.1 fig. The sharp peak given by 
the result of Eq. (^ will be broadened in an experimental 
situation. One of the main causes of broadening is that 
it is difficult to fix the particle number with high accu- 
racy. We have simulated this by assuming that the par- 
ticle number (fig) varies from experiment to experiment 
according to a Gaussian distribution (the width of the 
Gaussian is used as a measure of the fluctuations), and 
averaged over the different results. In these calculations 
we have introduced corresponding Gaussian fluctuations 
in the gap energy, because it is determined by the parti- 
cle number. For as large as 10% fluctuations in the gap 
the peak is still clearly visible although broadened and 
slightly shifted. 
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FIG. 2. The absolute value of the current / = —{Ne) as a 
function of the laser detuning A. The chemical potentials are 
fig — lOOkHz and fie = 0. The current is calculated assuming 
1%, 5%, and 10% fluctuations in the gap energy. Correspond- 
ing fluctuations in the particle number (fig) are taken into 
account using the exponential dependence between the Fermi 
energy and the gap energy. The unit of the current is arbi- 
trary. In the case of a normal state, there would be a sharp 
peak at A = 100 = fig — fXs- Due to the superconducting gap 
this peak is shifted and its shape modified. 



It is possible to prepare the desired initial number of 
atoms in states |e), \g) and \g') very accurately by apply- 
ing suitable incoherent laser pulses — this is one of the 
advantages of atomic Fermi gases compared to condensed 
matter systems. The number of atoms in the states \g) 
and \g') should be nearly equal to produce the supercon- 
ducting state, whereas the number of atoms in state |e) 
is not critical. In the simplest case the state |e) would 
be empty in the beginning. The use of our method then 
requires to know /ig rather accurately, since it appears in 
the tunable argument A = Afi + A = —fig + A of the cur- 
rent /(A). If this is difficult, it might be better to use an 
incoherent pulse to prepare initially Afj, = fi^, ^ fig — 0. 
In this case the chemical potential dependence disappears 
from the expressions in the homogeneous case. 

Finally, we discuss how to choose the states |e), \g) and 
\g') for a real atom (for instance ^Li). It will probably be 
difficult to find states in which the atoms do not interact 
at all. A more likely solution is to choose the states so 
that the interaction between all of them is weak and re- 
pulsive, except between \g) and \g') strong and attractive. 
Even a weak and attractive interaction instead of weak 
and repulsive would do, because the gap, critical temper- 
ature, and other essential BCS quantities have an expo- 
nential dependence on the scattering length, so the effect 
becomes easily negligible. Particularly ^Li seems to offer 
a variety of interaction strengths between different states. 
As shown in in a magnetic field of about B ~ 0.027 
the scattering lenght of several low field seeking hyperfine 
states is as strong as as ~ — 2000ao whereas the scatter- 
ing length between two high field seeking states is around 
zero. This is an example of the potential of modifying 
the scattering lengths by magnetic fields. In this case the 



trap should be optical in order to confine also the high 
field seekers. Alternatively, in a magnetic trap atoms in 
these states (corresponds to the non- interacting state |e) 
in our calculation) would simply fly out — this could even 
make the detection of the current easier. The detection 
of the atoms in the state |e), that is, the measurement 
of our basic observable /, could be done for instance by 
direct resonance fluorescence, or when a Raman process 
is used, indirectly from the intensity difference of the two 
Raman beams. 

In conclusion, we have proposed a method to observe 
the superconducting order parameter (gap) in a gas of 
cold Fermionic atoms. The idea is based on creating a 
normal phase - superconductor interface by coupling in- 
ternal states of different scattering lengths by a laser. 
The advantage of our scheme is the sensitivity to even 
small currents of atomic population caused by the laser 
interaction. Furthermore, our results extend beyond 
measuring the order parameter: they can be used as a 
starting point in describing various processes related to 
the transfer of atoms between two gas components (in su- 
perconducting or normal states) . If both of the internal 
states used were strongly interacting, one could perhaps 
observe Josepson tunneling of pairs from one supercon- 
ductor to another. With small modifications, our results 
also describe the functioning of an outcoupler from a su- 
perconducting gas of Fermionic atoms. 
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